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Abstract 

Experimental studies of FeCoB/MgO/FeCoB tunnel junctions indicate that boron diffuses into 
MgO during rf-sputtering and forms polycrystalline Mg-B-0 regions. These tunnel junctions pro- 
vide high tunneling magnetoresistance values and low RA products. However the crystal structure 
of the Mg-B-0 region remains unknown. Using density functional techniques, I examine three po- 
tential Mg(B) oxides including Mg2B20s (monoclinic and triclinic) and the orthorhombic mineral 
Kotoite (Mg3B20g). Kotoite is the best candidate for formation in magnetic tunnel junctions. The 
(100) surface of Kotoite has a good lattice match with (001) MgO and could template neighboring 
FeCo into bcc layers during annealing. Complex band structure analysis of Kotoite shows that 
the C2v Ai band has a much smaller imaginary k component than the C2 V A4 band. Based on 
symmetry analysis, the majority spin Ai band in FeCo should couple well with the Kotoite Ai 
band, while the minority FeCo A5 will couple partially with the A4 band. Kotoite provides a 
new route to high tunneling magnetoresistance based on spin filtering by a lower symmetry oxide 
region. 
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Magnetic tunnel junctions (MTJs) based on sputtered amorphous CoFeB leads and 
MgO tunnel barriers exhibit high tunneling magnetoresistance (TMR) at room tempera- 
ture (604%) [1J . High TMR values coupled with low RA products are essential for memory 
storage devices and CoFeB/MgO/CoFeB MTJs with low RA products have recently been 
demonstrated [2]. Boron ensures that deposited CoFe leads are amorphous, reducing in- 
terface roughness j3] and allowing MTJs to grow on synthetic antiferromagnetsptj. Upon 
annealing, MgO acts as a template for the growth of crystalline bcc FeCo layers on either 
side of the oxide [3] and provides the epitaxial interface required for spin filtering in the oxide 
layer [5J. 

While CoFeB leads have been very successful in improving MTJ properties, the exact 
role played by boron and where it resides after sputtering and subsequent annealing has 
remained a mystery. Originally it was assumed that all boron diffused away from the MgO 
layer during annealing. This would leave isolated FeCo layers at the MgO interface, resulting 
in high TMR values. However, recent work indicates that the movement of boron can depend 
strongly on deposition and annealing conditions. A theoretical study found that a single 
layer of boron at the CoFe/MgO interface was energetically favorable and could significantly 
reduce TMR[7j. However, this work did not examine the possibility of boron diffusion and 
incorporation into MgO. Recent experimental studies provide evidence that boron oxide 
does form during sputtering deposition [HI [91 [10] and that a significant amount of boron can 
diffuse into MgO during annealing and enhance TMR for thin oxide tunnel junctions[llJ. 
Mg and O K-edge EELS data from these studies indicated that Mg and O coordination was 
different from that of bulk MgO, indicating the presence of Mg-B oxides. Boron K-edges in 
the oxide region show that boron is oxidized with a BO3 configuration that is present after 
deposition and remains even after annealing [12]. 

Experimental studies of MTJs with Mg-B-0 oxide regions show that boron incorporation 
into the oxide region (1) preserves a high TMR value upon annealing, (2) provides MTJs 
with the low RA product needed for industrial applications, and (3) acts to template the 
neighboring FeCo leads to grow in bcc crystalline layers. However, a clear understanding of 
the specific Mg-B-0 crystal structure in this region is crucial for understanding the under- 
lying physical properties of this system and how they may be optimized. I examine how the 
incorporation of boron affects the local density of states and use complex band structure 
analysis to address whether symmetry based spin filtering is still possible in these Mg-B-0 
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oxides. 

In this Letter, I investigate potential boron based oxide formation within tunnel junc- 
tions. I consider three candidate Mg-B oxides including magnesium borate (Mg 2 B 2 05) in its 
monoclinic [13] (space group P2i/c) and triclinic[H] (space group PI) forms as well as the 
mineral Kotoite (Mg 3 B 2 06) (space group Pnmn)P2] (see Fig. [TJ. In each crystal, the Mg 
atoms are surrounded by an octahedron of oxygen atoms, shown in green in Fig. [1} Boron 
atoms are always surrounded by three oxygens, BO3, in a triangular configuration shown as 
blue triangles. The main difference between Mg 2 B 2 0s and Mg 2 B 3 06 comes from the coor- 
dination of BO3 triangles. In Kotoite (Mg 3 B 2 0e), single BO3 triangles link chains of oxygen 
octahedra, while in Mg 2 B 2 0s crystal structures, double BO3 triangles link together oxygen 
octahedra[16j. While the monoclinic form of Mg 2 B 2 0s is found in nature as the rare mineral 
Suanite[13], synthetically produced magnesium borate is used in industry as a catalyst for 
hydrocarbon conversion [17] and for friction reduction [TS]. Monoclinic and triclinic Mg 2 B 2 05 
nanowires have also been grown |2*D] . 

I examined the impact of boron on the density of states (DOS) in the oxide region. Density 
functional calculations were performed using both a pseudopotential-plane wave approach 
(as implemented in Quantum Espresso[2T]) and a full potential augmented spherical wave 
approach [22 j. Since very little work has been done on the electronic structure of these sys- 
tems, using multiple ab-initio techniques provides an additional accuracy check. All calcu- 
lations were done within the local density approximation. A 100 Ryd cutoff energy was used 
for plane wave calculations. Monkhorst-Pack k-point grids of 12x12x12, 8x16x8, 10x6x12, 
and 8x4x12 were used for MgO, Mg 2 B 2 0s (monoclinic), Mg 2 B 2 0s (triclinic), and Kotoite, re- 
spectively. Mg and B ions were described using von Barth Car pseudopotentials|23j. while an 
ultrasoft pseudopotential was used for O atoms [21]. Exchange and correlation interactions 
were represented by the Perdew and Zunger parameterization |25j . Lattice constants and 
atomic positions were relaxed based on Hellman-Feynman forces. For augmented spherical 
wave calculations, empty spheres were added to insure a proper description of open regions 
in the crystal structure. 

Introducing boron into the MgO region will affect the density of states in the valence 
and conduction band. In Figure [2j the DOS of the three Mg-B-0 systems is compared to 
MgO. The MgO conduction band DOS gradually decreases as it approaches the conduction 
band edge. In contrast, all three Mg-B-0 systems show two sharp peaks at the conduction 
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(c) (d) 

FIG. 1: (color online) The crystal structures for (a) MgO, (b) Mg 2 B 2 05 (monoclinic), (c) Kotoite 
Mg3E>206, and (d) Mg 2 B 2 05 (triclinic). Oxygen atoms are given in red and magnesium atoms are 
located at the center of the green octahedra. Boron atoms are located at the center of the blue 
triangles. 

band edge. In addition, the valence s and p band DOS in the monoclinic and triclinic forms 
of Mg 2 B 2 05 are split with the formation of two sharp secondary peaks at slightly lower 
energies denoted by arrows. While the valence DOS for the s and p bands in Kotoite is 
slightly broadened in comparison to MgO, separate DOS peaks do not form. 

The cause of these DOS changes for MgB oxides can be determined by examining the 
partial DOS resolved into angular momentum channels (Fig. [3]). The sharp peaks at the 
conduction band edge in all three Mg-B-0 crystal structures are due to boron p states. In 
monoclinic and triclinic Mg 2 B 2 05, the splitting of the s and p valence DOS is due to the 
formation of sharp peaks related to boron-oxygen bonding. The monoclinic and triclinic 
forms of Mg 2 B 2 05 are characterized by double boron-oxygen triangles which are connected 
by an oxygen in the center. This oxygen has a different bonding environment than the 
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FIG. 2: (color online) The total density of states is shown for (a) MgO, (b) M^E^Os (monoclinic), 
(c) Kotoite Mg3B2C>6, and (d) Mg^E^Os (triclinic). The Fermi energy level is represented by a 
vertical dashed red line. The additional valence band peaks which form in the monoclinic and 
triclinic forms of Mg2E>205 are denoted by blue arrows. 



other oxygens that form the Mg octahedras. This leads to splitting of the valence band 
DOS. Kotoite, in contrast, only has a single BO3 triangle and lacks this DOS splitting in 
the valence bands. 

Although EELS data shows boron forming a relatively uniform distribution of BO3 tri- 
angles in the MgO region [12], experimentally it is difficult to resolve whether this leads to 
Mg-B-0 crystals or clusters of B2O3 in a MgO matrix. To rule out B2O3 cluster formation, 
I also compared the total energy of the different Mg-B-0 compounds to the formation of 
separate MgO and B2O3 regions. The trigonal crystal structure for B203[26] was used for 
total energy calculations. The energy difference per atom in comparison to separate B2O3 
and MgO regions is -0.0944 eV for Kotoite, -0.1005 eV for monoclinic Mg 2 B 2 5 , and -0.1006 
eV for triclinic Mg 2 B 2 05. In all cases, the Mg-B-0 crystal structures are more stable than 
separate regions of MgO and B2O3. This indicates the formation of isolated B2O3 clusters 
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FIG. 3: (color online) The partial density of states is shown for (a) monoclinic Mg2B20s and (b) 
Kotoite Mg2B3C>6. The partial density of states is only plotted for inequivalent atoms. The Fermi 
energy level is denoted by a red dash line. 

is unlikely and that the oxide region is made up of a crystalline magnesium boron oxide. 
However, all three structures are fairly close in terms of total energy and other factors must 
be considered to determine the correct Mg-B-0 structure. 

Several key pieces of information support Kotoite formation in rf-sputtered MTJs. The 
measured boron concentration in these Mg-B-0 regions is fairly low (5-10%) [HI This 
favors a Mg-B-0 crystal with a low boron concentration. Of the three systems considered 
in this work, Kotoite, with single BO3 triangles, has the lowest boron concentration (18%). 
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Experimental studies also indicate that the Mg-B-0 region is crystalline after annealing and 
has a good epitaxial match with MgO on one side and still acts to template the neighboring 
FeBCo leads into bcc crystalline layers. This indicates that the Mg-B-0 crystal has a good 



lattice match with MgO in the (100) direction. In Figure 4(a) , the MgO (100) crystal face is 



shown with Mg atoms in green. The (100) face of Kotoite is shown in Figure 4(b) Kotoite 
is orthorhombic and the (100) face is rectangular. The lattice match with (001) MgO in 
this direction is very good. The B lattice vector of Kotoite is 8.416 A which is very close to 
twice the MgO lattice constant (2Am 9 o — 8.422A). The C lattice vector of Kotoite is 4.497 
A which is slightly larger than the MgO lattice vector, 4.211 A. The Mg atomic positions 
on the (100) face of Kotoite also mimic those found in (001) MgO for almost all sites. The 
monoclinic and triclinic forms of Mg2B20s in contrast provide a poor lattice match to MgO 
and would be unable to template FeCo into bcc layers. 

Zhang and Butler et al. showed that the high TMR in FeCo/MgO/FeCo MTJs relies 
on symmetry selective filtering in MgO|2H]- In particular, the FeCo majority carrier A x 
band couples efficiently with the Ai complex band in MgO, while minority carriers can only 
couple through the A5 and A2 bands which decay much faster. In order for Kotoite to provide 
high TMR values, it must provide similar symmetry selective filtering to that of MgO. To 
determine whether this is the case, the complex band structures for both the (001) direction 
in MgO and the (100) direction in Kotoite were calculated using the approach of Choi and 
Ihm[2"9"] as implemented in the Quantum Espresso package [3Q]. The complex band structure 



for the tetragonal unit cell of (001) MgO is shown in Figure 4(a) In the (001) direction 
for MgO, the band symmetry group is C± v . As has been reported in previous studies, there 
is a substantial difference in the imaginary component of the Ai and A5 complex bands 
that leads to efficient majority spin filtering. The complex band structure for Kotoite in 
the (100) direction is shown in Figure 4(b)[ The band symmetry group along the T to X 
line is Ci v . Two different imaginary bands labeled Ai and A 4 intersect the Fermi energy 
at -0.28 and -0.46 respectively. These imaginary bands are continuations of real conduction 
bands based on boron p states. The imaginary component of the two bands determines 
how effectively electrons tunnel through the barrier. The large separation in the imaginary 
components in these two bands indicates that electrons that can couple to the orthorhombic 
Ai band should dominate transport through the oxide region. Introducing FeCo leads on 
either side of the oxide region will cause the position of the Fermi energy level to change. 
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FIG. 4: The complex band structure of MgO and Kotoite (MgsE^Og) are shown in panels a and b, 
respectively. K-vector components are shown in terms of positive and negative regions of the plot 
to indicate real and imaginary components, respectively. Real bands are denoted by black lines. 
Purely imaginary bands are given by blue lines and complex bands are denoted by orange lines. 
The (001) MgO surface and the (100) Kotoite surface are shown as insets in panels (a) and (b), 
respectively. Mg atoms are denoted by green circles. 

However, over the entire energy range in the band gap, the imaginary components of Ai is 
significantly smaller than A 4 . 

This interface presents an interesting problem in terms of symmetry filtering for MTJs. 
The FeCo bcc layers and (001) MgO all possess cubic symmetry [6]. However, Kotoite is 
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orthorhombic and possesses C2 V symmetry along the [100] direction. This raises the question 
of how the cubic symmetry states (A 1; A 5 , and A 2 ) will couple into the A x and A 4 bands 
in this lower symmetry region. The cubic C4„ A x bands transform as linear combinations 
of functions with 1, z, and 2z 2 — x 2 — x 2 . The Ai bands in orthorhombic systems 
transform as linear combinations of 1, z, x 2 , y 2 , and z 2 . Of the symmetry states available in 
Kotoite, the Ai band should couple most effectively with the majority spin Ai band in the 
FeCo magnetic leads. For the minority carriers, the cubic A 5 band transforms as linear 
combinations of zx and zy symmetry. The closest analogue in the orthorhombic system is 
the A 4 band which transforms as yz and has a much larger imaginary component than the 
orthorhombic Ai band. Based on this symmetry analysis and the complex band structure, 
Kotoite should serve as an effective spin filter similar to MgO. 

In this Letter, I have used density functional approaches to examine three possible magne- 
sium boron oxides that could form during rf-sputtering of magnetic tunnel junctions. Kotoite 
(Mg 3 B 2 06) is shown to be a strong candidate for symmetry based spin filtering. Density 
of states analysis shows that all three compounds have sharp boron p bands situated at 
the conduction band edge. Several factors indicate that Kotoite is the best candidate for 
formation in MTJs. The boron concentration in this structure is closest to that measured 
in experiments. Kotoite grown in the (100) direction also provides a good epitaxial match 
to MgO (001) and can act to template neighboring FeCo regions as crystalline bcc layers. 
The diffusion of boron into the MgO region may also help passivate Mg or O vacancies 
that would otherwise reduce the TMR in the magnetic tunnel junctions[2D |32]. Complex 
band structure analysis shows that the cubic (001) Ai should couple with the orthorhombic 
(100) Ai bands in Kotoite, while the cubic (001) A5 bands should partially couple with 
the orthorhombic (100) A4 bands. There is sufficient difference in the decay rates of the 
imaginary bands so that spin filtering should be expected. This Letter has focused on the 
potential formation of Kotoite under present MTJ fabrication conditions. Further experi- 
mentally research is required to determine conclusively whether Kotoite is indeed present in 
the barrier region. Regardless of this issue, the current study shows that an experimental 
effort focused on forming MTJs based on Kotoite barriers is worthy of exploration. Using a 
lower symmetry oxide region between two cubic symmetry magnetic leads could provide a 
new route for enhanced spin filtering and potentially a new generation of high TMR devices. 

Calculations were done on the Intel Cluster at the Cornell Nanoscale Facility, part of the 
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